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Abstract 

We confront two types of phantom dark energy potential with observational data. The models 
we consider are the power-law potential, V oc (fi^, and the exponential potential, V oc exp (Xcft/Mp). 
We fit the models to the latest observations from SN-Ia, CMB and BAO, and obtain tight con- 
straints on parameter spaces. Furthermore, we apply the goodness-of-fit and the information 
criteria to compare the fitting results from phantom models with that from the cosmological con- 
stant and the quintessence models presented in our previous work. The results show that the 
cosmological constant is statistically most preferred, while the phantom dark energy fits slightly 
better than the quintessence does. 
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I. INTORDUCTION 



Observations over the past dozen years have shown that the universe is currently under 
accelerating expansion (see for reviews). Under the framework of general relativity 

and standard cosmology, a new form of exotic energy with negative pressure (p < — p/3) is 
required to explain this phenomenon. Current observations suggested this so called "dark 
energy" made up about 73% of the energy density of the universe js] 0| 0. So far wde 
has been constrained to be very close to —1 assuming it is constant and the universe is 
flat. This seems to suggest the observation data prefer a cosmological constant as dark 
energy. However, a parametrized dark energy w = w + w a (l — a) has ~ 100% uncertainty 
in dynamical parameter w a , when compared to data. This means dynamical dark energy 
models are not excluded. 

Several dark energy model have been proposed in order to explain the cosmic acceleration. 
In addition to the most discussed cosmological constant, dynamically evolving scalar-field 
dark energy has been widely studied (for examples, see p| [3| [8j). Quintessence is a specific 
case of scalar- field dark energy with canonical kinetic terms, which admits — 1 < < 1. 
It has drawn much attention because it can in principle provide the "tracker" property - 
a property that makes the energy density today insensitive to its initial condition, i.e., the 
initial energy density value js], (9]. While fine-tuning of potential parameter is still 



required 10], tracker quintessence can alleviate the cosmic coincidence problem because a 
precise setting of initial energy density ratio for matter and dark energy is no longer required. 

Present observational constraint on w^g still allows for %g < — 1. For example, results 
from WiggleZ [5] showed current constraints on constant w is —1.114 < w < —0.954 after 
combining with the latest SN-Ia, CMB and BAO data. We note that while quintessence 
only allows for— 1 < < L, phantom scalar-field, another dynamical DE model, proposed 
by Caldwell, Carroll et al Uj, that invokes a negative kinetic energy, can satisfy < — 1. 
While there exist several known theoretical difficulties for phantom scalar-field dark energy 
model such as the violation of null dominant energy condition (NDEC) 1 1 21 and a possible 

nn 

Big Rip phase in the future [13J] [14J, it is still very worth while to confront it directly and 
independently against the cosmological observations, especially since the current best-fit for 
equation-of-state is smaller than — l(see JiJ [l^] 16] for examples). 

In our previous work [17|, we have tested several tracker quintessence models with ob- 



servational data. The result showed that the best-fit of the inverse-power-law potential and 
the inverse-exponential potential models both reduced to the cosmo logical constant. Moti- 
vated by this and the implication of wde < ~~ 1 from observations, in this paper we put the 
phantom dark energy models to the test. We consider two specific scalar- field phantom po- 



tentials: the power-law potential 



15] and the exponential potential 18] . each with one free 



parameter. The reason to choose these potentials is that they also possess the attractor-like 
property: insensitive to initial conditions. We take the model-based approach to confront 
the models with observational data. The data we use includes the latest Type-la supernova 
(SN-Ia) compilation set, the cosmic microwave background (CMB) and the baryon acous- 
tic oscillations (BAO) observations. We also confront these models with the cosmological 
constant and the quintessence scalar-field models by using the Goodness-of-Fit test and the 
information criteria to assess the merit of each model. 



II. TRACKER AND ATTRACTOR PHANTOM 
A. Phantom Formalism 

Phantom dark energy with equation of state w < — 1 is achieved by introducing a negative 
kinetic energy term in the action. In this way phantom scalar field is slowly "rolling up" 
the potential. The energy density and pressure of phantom field can then be given as 

P^ = ~<P 2 + V(<P), (1) 

P* = -\?-V(<l>). (2) 

with w<f,(z) = p<j>/P(f> < —1 in the range < 2 /2 < V(<f>). The evolution of the phantom 
field is governed by its equation of motion: 

in which H denotes the Hubble expansion rate, the dot denotes the derivative w.r.t. the 
physical time. Assuming a flat universe, the Friedmann equation can be written as 

H 2 (z) = (p r {z) + p m (z) + p+(z)) 



Hi 



Q r {l + zf + n m (l + zf + % exp ^3 J [1 + w^z')] 



(4) 
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where p r (z) is the radiation energy density, p m {z) is the matter energy density, Ptj>(z) is the 
scalar field energy density, H is the Hubble constant. 



B. Tracker Phantom 

In the quintessence scenario, there exists a special "tracker solution" to which other 
solutions would converge j^] . A wide range of initial conditions for and <p will approach 
a common evolutionary track of p^ and w^; this means the tracker field model is insensitive 
to its initial conditions. A very large range of initial values of p\ is thus allowed without 
changing cosmic history. This property makes the tracker model very interesting to study, 
because it can alleviate cosmic coincidence problem. Conditions for tracker quintessence are 
such that T = V"V/(V') 2 > 1 and is nearly constant (\V /T(V /V)\ < 1) for a wide range 
of plausible field initial conditions. Here the prime denotes the derivative w.r.t <fi. Under 
these conditions, the tracker solution for quintessence can be approximated as 9| 



w *~ i + 2(r-i) ' (5) 



where wb is the background dominant component in the equation-of-state. The fact that 
T > 1 ensures that < wb, so that at late times dark energy will eventually take over and 
become dominant. 



The tracker solution for phantom has been studied in 19j 20[. Its tracker condition 
is closely related to that for quintessence: T = V"V/(V') 2 < 1/2, and is nearly constant. 
Because < — 1 all the time, the energy density of phantom dark energy either stays the 
same or grows with time. This ensures that the dark energy density will eventually take 
over and become the dominant substance at late times. The form of its tracker solution is 
the same as in Eq. (jSJ). It is evident that the tracker phantom models are also insensitive to 
initial conditions for phantom scalar field. One example of tracker phantom is a power-law 
potential V oc M with < p < 2. In Fig. (TJ we show an example for the tracking behavior 
of the power-law potential phantom. 




FIG. 1. Example of tracker phantom, with power-law potential. 
C. Attractor Phantom 



The tracker solution is just one type of "attractor solutions" for scalar fields 19] 21 ]. 
Attractor solutions are the stable critical points of autonomous equations (re-written from 
the equation of motion and Einstein equation) to which different initial conditions converge 



(see proof in 19] 21( 22j). In contrast, the tracker solution is not a usual attractor because 
its critical point is not fixed, but changes with time when the background fluid dominates. 
Other than the tracker, there are two additional types of late time attractor solution for the 
phantom scalar field: the Big Rip attractor and the de Sitter attractor. These two attractors 
can only be reached in the future (when = 1); however we found numerically that the 
exponential potential model in the Big Rip attractor case remains insensitive to some range 
of the initial values of <fi and 0, as long as fa <C Mp. This allows certain range of initial p^, 
and is thus still worth looking into. 

The Big Rip attractor acquired its name because the attractor solution approaches < 
— 1 in the future. This will cause a catastrophic "big rip" , where the energy density and 



the scale factor will blow up within a finite time 13| 23[. Although this may seem an 



unappealing feature, it is theoretically permissible and is therefore worthy of investigation. 
The condition for the Big Rip attractor is T ~ 1 and A = —V'/kV ^ [19J], k 2 = 8nG. One 
example is the exponential potential V oc exp (0/Mp) mentioned above. 

In this paper, we consider two phantom scalar-field dark energy models: the power-law 
potential V oc </> M and the exponential potential V oc exp(A0/Mp). For the power-law 



potential phantom, it has the tracker property. For the exponential potential phantom, 
there exists an attractor solution that approaches the big rip in the future. We analyze 
these two models with observational data, obtain constraints on the model parameter, and 
assess the merit of the models derived from the best-fit results. 



III. DATA 



We use observational data from SN-Ia, CMB and BAO described below. 



A. Type-la Supernovae 



We use the latest SNIa dataset, Union2.1 compilations, released by the Supernova Cos- 
mology Project that contains 580 SN-Ia in the the redshift range of 0.02 < z < 1.5 j^]. 



24| |32|. The dataset provides the distance 



This compilation includes supernova data from 
modulus that contains information of the luminosity distance that can be used to constrain 
the dark energy. 

The distance modulus is defined as following: 



H th (z) = 5 log 



10 



d L {z) 
Mpc 



+ 25 = 51og 10 (D L (z))+fi , 



(6) 



where Dl(z) = Hodi(z) is the Hubble-free luminosity distance. We marginalize XsNia over 
the nuisance parameter /i by minimizing it with respect to /io- The marginalized XsNia 



is 



33|-|35| 



where 



XsNia — A "^7> 



B 2 



(7) 



A = ^2 [ 5 D L(zi, par) - n ohs (zi) C^ 1 5 log 10 D L (z j} par) - /j, obs (zj 
B = [ 5lo &io D dzi,par) - /i obs {zi) C£\ 



(8) 

(9) 
(10) 
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B. Cosmic Microwave Background 



The 7- year WMAP data provides the "distance prior" that can be used to constrain dark 



energy (4]. The distance prior includes the CMB shift parameter R = 1.725 ± 0.018 given 
by 



R= ^n m H$(l + z*)D A (z*), (11) 
and "acoustic scale" I a = 302.09 ± 0.76 given by 

U = (l + z.) A ( y *\ 12 

where z* is the redshift at decoupling, Da is the physical angular diameter distance, and r„ i s 



the comoving sound horizon. We use the fitting formula proposed by Hu and Sugiyama 36] : 



z* = 1048 [1 + 0.0012A(Q b h 2 )- - 738 ] [l + gi {Q m h 2 ) 92 } , (13) 

_ 0.0783(n^ 2 )- - 238 

91 ~ l + 39.5(fi 6 /i 2 ) - 763 ' ^ ' 

0-560 

92 l + 21.1(JVi 2 ) 1 ' 81 ' 

The comoving sound horizon is 

/ \ 1 f 1/(1+Z) da 
r (z) = / ^^^^^^^^= (16) 

V3Jo a?H{a)^l + {m b /m i )a 

where Q b is the baryon density and f2 7 is the photon density. We construct Xcwib = 
^2ij[ x i ~ x ? bs ]^'ij 1 [ x j ~ x f bs ], where C^ 1 is the inverse covariance matrix given in |4j, and 



C. Baryon Acoustic Oscillations 

We followed [5| and use three sets of BAO distance dataset: 6dFGS, SDSS and WiggleZ, 
for our study. 

We use the joint analysis of the Two Degree Field Galaxy Redshift Survey (2dFGRS) 



data [371] and the Sloan Digital Sky Survey (SDSS) Data Release 7, which provides two 



distance measures of d .35 = r s (z d )/ D v (0.35) = 0.1097±0.0036 and g? . 2 = r s(z<i) I 'Dy (0.2) 
, where r s (z d ) is the acoustic sound horizon at the drag epoch, D v 



0.1905 ± 0.0061 



1 /3 

[(1 + z) 2 D 2 A (z) I H (z)] . The fitting formula for Zd is defined by Eisenstein & Hu [39j. The 
Xbao! is I2ij[ d i ~ df^C^ldj ~ d f% where d i = (^0.2,^0.35) and 



c- 1 



30124 -17227 
-17227 86977 



(17) 



The fitting formula for Zd has the form: 

1291(ft m /i 2 ) - 251 



1 + 0.659(fi m /i 2 ) - 828 



(18) 



where 



61 = 0.313(fi m /i 2 )- a419 [1 + 0.607(fi m /i 2 ) - 674 ] , 6 2 = 0.238(fi m /i 2 ) a223 . (19) 

The second is the 6dFGS data in [40j, which provides 4(0.106) = 0.336 ± 0.015. We thus 
have x 2 bao 2 = [4(0.106, par) - 0.336] 2 /(0.015) 2 . 

Finally, we include the result from the WiggleZ Survey |5[ . WiggleZ provides three corre- 
lated measurements: A obs = (A(z = 0M),A(z = 0.6), A(z = 0.73)) = (0.474,0.442,0.424), 
with the inverse covariant matrix 



/ 



c- 1 



1040.3 -807.5 336.8 
3720.3 -1551.9 
2914.9 



\ 



(20) 



and A(z) defined as 



A{z) 



cz 



(21) 



The Xbao 3 can ^ e wrrfc ten as [A° bs — A^C^lAf 8 — Aj]. It is obvious that Xbao = Xbao 1 + 



2 2 

Xbao 2 + Xbao 3 



D. Prior 



For the radiation, we fix Q 1 = 2.469 x 10~ 5 /h 2 , and use the relation Q r = + 
0.2271N e ff) as the radiation e nerg y density 41J. N e ff is the effective number of neutrino 
species and is taken to be 3.04 [41]. We further impose the prior of H = 73.8 ± 2.4 km s _1 
Mpc" 1 from {4^. This prior is an independent and complementary constraint on parameter 
h. The total chi-square Xtotai = XsNia + Xcmb + Xbao + Xh is marginalized overtt b h 2 and 



the reduced Hubble constant h by minimizing Xtotai with respect to Qbh 2 and h [41]. 



IV. OBSERVATIONAL CONSTRAINTS ON TRACKER AND ATTRACTOR 
PHANTOMS 



We consider two potential forms, the power-law potential V — V\ (</>/Mp) M and the ex- 
ponential potential V = V2 exp (A0/Mp). As mentioned above, the power-law potential 
corresponds to the tracker phantom, whereas the exponential potential admits a late time 
big rip attractor solution. Here [i and A are positive, dimensionless constants. The constant 
V\ and V2 are determined by requiring total energy density today equals to the critical en- 
ergy density p c in a flat universe (i.e Q to t — !)■ We calculate our x 2 by solving Eq. fl3]) and 
Friedmann equation Eq. (jlj) numerically. 

The results are given in Table [B For comparison, we also provide results of two quintessence 
models, the inverse-power-law potential and the inverse-exponential potential. The best- 
fit equation of states for two phantom models at late times are given in Fig. |2l and the 
two-parameter likelihood contours ((fi m , a), (Q m , A)) are given in Fig. [3j 

In order to test the merit of the model, we perform the goodness-of-fit (GoF) test to all 
models. The meaning for GoF is, assuming a model to be true, the probability of finding a 
new set of data that gives worse \ 2 than that deduced by the current data. The higher the 
GoF is, the more viable is the model. Explicitly, it is defined as T(i//2, x 2 /2)/r(z//2), where 
r(V/2, x 2 /2) is the upper incomplete gamma function and v is the degrees of freedom. 

To further assess the relative merit between models, we invoke the "information crite- 
ria" . The information criteria (IC) is a set of statistical considerations that take both data 
fitting and model complexity into account; they favor models with better fit and fewer pa- 
rameters. The application of the information criteria to cosmology was first launched by 



Liddle 43]. Here we consider two kinds of information criteria: the Bayesian information 



criterion (BIC) 44] and the Akaike information criterion (AIC) [451 ] . The BIC is given by 
BIC = —2 In C max + k In N, where C max is the maximum likelihood, which is equivalent to 
the minimum x 2 for gaussian errors, k is the number of parameters, and N is the number 
of data points used in the fit. The second term klnN serves as the "penalty" to the model 
that invokes extra parameters. The AIC is defined as AIC = —2 In C max + Ik. Both BIC 
and AIC favor smaller values while BIC charges stiffer penalty for extra parameters when 
In TV > 2. Taking the BIC and AIC for the cosmological constant as the reference value, 
we compute the differences in BIC (ABIC) and AIC (AAIC) of other dark energy models. 



They are listed in Table [I] 



TABLE I. Fitting Results 



Model 


Best-fit x 2 


Best-fit parameters^ 


GoF ABIC AAIC 


Cosmological Constant 


Xlt = 550.41 

XSNIa = 545 - 24 

XImb = 0-19 
Xbao = 2-32 


O _n9S1+ - 013 
i4 m — <J.z,o±_ 0.012 


85.8% 0.0 0.0 


Phantom with 

V oc exp (Xcp/Mp) 


XL = 550.09 
xk/„ = 545.66 
XImb = 0-37 
XIao = 2-06 


O n OQfi+0.013 

ic "m — u,ZlOU — 0.012 

A = 2.99+ 210 


85.4% 6.1 1.6 


Phantom with 

V OC (fr 


xL = 550.17 

XSNIa = 545 - 53 

A ,2 n a 1 

XcA/B - U - 4i 

Xbao = 2-17 


o n oan+o.oi3 

ai m — u ou — 0.012 

,, — n nsi +0.154 

/i — U.Uol 


85.3% 6.1 1.8 


quintessence with 

V oc (j)- a 


Same as cosmological constant 


O _n9S1+ - 013 
iL m — u.z,ox_ 0.012 

a = 0+ - 11 


85.1% 6.4 2.0 


quintessence with 
V oc exp(/3M P /0) 


Same as cosmological constant 


O _n981+ - 013 
iL m — u.z,ox_ 0.012 

/3 = 0+ 0010 


85.1% 6.4 2.0 



the best-fit and the 68.3% confidence interval for each parameter (A% 2 = 1) 



V. DISCUSSION 

We have tested two potential forms of phantom scalar-field dark energy, the power- 
law V = V\ {^jMpf and the exponential potential, V = V 2 exp {Xcp/Mp), with current 
observations. Tight model parameter constraints are obtained in Table [B We also provided 
results of cosmological constant and two types of quintessence potentials, the inverse-power 
law potential and the inverse-exponential potential. Cosmological constant yields the best 
goodness-of-fit and smallest information criteria. This shows that the cosmological constant 
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-0.96 - 



■0.98 - 



-1.02- 



-1.06 



V oc exp(?i4»/M p ) 
Cosmological Constant 



J I I L 



J I I l_ 



0.5 



1 
Z 



1.5 



--0.96 



— 0.98 




j i i i_ 



■1.06 



FIG. 2. The late time evolution of w^iz) corresponding to the best-fit parameters. Solid curve 
corresponds to V oc (ff 1 , dot-dash curve corresponds to V oc exp (X(f>/Mp). 



is still the most preferred dark energy model among all that we have considered, even with 
the constraint w > — 1 removed. Phantom models fit worse to observations than that with 
the cosmological constant, but are slightly better than the two quintessence potential models. 
Although the current observations still prefer the cosmological constant as the dark energy, 
phantom and quintessence models under considerations are only slightly worse in terms of 
GoF and the information criteria; all our models in Tabled] have GoF~ 85%, ABIC~ 6, and 
AAIC- 2. 

Another interesting result is the best-fit for phantoms. Unlike quintessence, which have 
the best-fit equivalent to ACDM jl7], that for phantom models moves away from — — 1 
(see Fig. [2]), as the best-fit for parameters A and /x are no longer zero (Fig. [3]). For V = 

V x (0/Mpf , the best-fit w^O) 1.03; as for V = V 2 exp (\</>/M P ), we find w^(0) ~ -1.06 

for the best fit. This indicates that the current observations may prefer woe < — 1- 
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FIG. 3. Joint constraints on (fi m ,/i) and (f2 m ,A). The dark gray and the light gray regions 
correspond to the 68.3% and 95.4% confidence regions, respectively. 



In summary, the model-based approach we used in this paper suggests that the cosmo- 
logical constant is more preferred, and dark energy models with w < — 1 is preferred over 
w > — 1. Notice that in both phantom models the cosmological constant cases are still inside 
1 — er range. This means currently we cannot distinguish small-dynamical dark energy from 
the cosmological constant. Future observations from next generation dark energy probes are 
expected to constrain w about ten times better than the present value 46] . More stringent 
constraints on the parameter space are thus expected to be obtained. By then we should be 
able to attain more insights into the physics of dark energy models with this model-based 
approach, or even rule out some of the models at a sufficient confidence level (see results 
with projected data in 47] 48] (43 1 (sQ|). 
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